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Abstract 

Time evolution of a nonequilibrium polariton condensate has been studied in the 

■ framework of a microscopic approach. It has been shown that due to polariton-polariton 

■ scattering a significant condensate depletion takes place in a comparatively short time inter- 

& ; 

val. The condensate decay occurs in the form of multiple echo signals. Distribution-function 
dynamics of noncondensate polaritons have been investigated. It has been shown that at the 

o : 

O ■ initial stage of evolution the distribution function has the form of a bell. Then oscillations 

> : 

arise in the contour of the distribution function, which further transform into small chaotic 
ripples. The appearance of a short-wavelength wing of the distribution function has been 
demonstrated. We have pointed out the enhancement and then partial extinction of the 
sharp extra peak arising within the time interval characterized by small values of polariton 
condensate density and its relatively slow changes. 
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1 Introduction 



In a great number of works concerning spatial and time evolution of laser radiation in resonant 
media the classical description is used when the photon and material fields are considered as 
interacting macroscopically occupied coherent modes (condensates). When a quasiparticle quits a 
condensate mode, that is usually taken into account by introduction of phenomenological damping 
constants. In this case it is assumed that (1) the condensate decay is exponential and (2) the 
condensate time life is considerably longer than the duration of the coherent nonlinear process 
being investigated. 

From our point of view it is more consistent to use a microscopic approach (see, e.g., [[[], §. ||). 
where the dynamics of interacting condensate modes and quantum fluctuations appearing in a 
system as a result of real scattering processes of quasiparticles are described on equal footing. In 
other words, one should take into account (1) the interaction between condensates, (2) excitation 
of noncondensate particles, and (3) the backward influence of non-condensate particles on the 
condensates. 

The microscopic approach is used in this work to describe the of time evolution of the system 
of initially coherent high-density excitons and photons. 

It is known 0] that coherent electromagnetic radiation resonant to the isolated exciton en- 
ergy level excites in the crystal the coherent polariton wave with the wave vector ko 7^ — the 
nonequilibrium polariton condensate. Different scattering processes accompanying its propaga- 
tion lead to the loss of initial coherence of the polariton wave, complete or partial depletion of 
condensate, excitation of polaritons characterized by some statistics with wave vector k 7^ k , and 
other phenomena. 

At sufficiently high excitation energies the processes of polariton scattering due to exciton- 
exciton collisions prevail. This scattering mechanism is of considerable interest due to the recent 
experimental investigations |J and many interesting physical results obtained in theoretical study 
of dynamic and kinetic processes in the system of interacting polaritons. 

According to || [7j exciton-exciton scattering is of significant importance in the situation 
where coherent polaritons are excited in a certain spectral region where energy and momentum 
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conservation laws allow real processes of two-quantum excitation of polaritons from the condensate. 
These processes lead to instability of the completely condensed state of the polariton system. 
Existence of this spectral region situated around the isolated exciton resonance is due to the 
peculiarities of the polariton dispersion law. 

In H [7J the energy spectrum of non-condensate polaritons, arising as a result of decay of 
coherent polariton wave, is studied. According to || |7J in some regions of k space the energy 
spectrum does not exist. 

It should be mentioned that investigations performed in || |7j are based on the model for- 
mally analogous to that used by Bogoliubov in |§ to study the equilibrium system of weakly 
nonideal Bose gas. In the nonequilibrium situation considered in || [7|, when decay of the polariton 
condensate and excitation of noncondensate polaritons take place, this model is adequate to the 
real situation only at the initial stage of the condensate decay when the number of polaritons in 
the condensate is still much greater than the total number of non-condensate polaritons. But this 
stage is essentially nonstationary and the occurrence of condensate instability provides evidence 
for that. The study of the energy spectrum implies the determination of the steady states of the 
system || . That is why the results of [§, [7j , concerning the energy spectrum of the system, based 
on the above mentioned model, in our opinion cannot be regarded as well grounded. This remark 



refers to the relatively recent works [10|, [IjJ as well. 

Because of essential nonstationarity of the processes in the system methods of nonequilibrium 
mechanics should be used to describe it adequately. Derivation of the equations that describe 
kinetics of the polariton condensate decay and excitation of quantum fluctuations has some specific 
features owing to degeneracy in the system. As the energy and resulting wave vector of two non- 
condensate polaritons can be equal to the energy and wave vector of two condensate polaritons, 
respectively, there is degeneracy of two-particle states. Moreover, the presence of the condensate 
in the system also leads to degeneracy due to its macroscopic amplitude fL2| . 

The correct description of the system with degeneracy requires the introduction of abnormal 
distribution functions [HJ 

*k = <*ko> 
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and 



F k = <$ k $2ko-k> 



together with the normal (usual) ones 



N k = ($+$ k ). 



Here $j^($k) ar e Bose operators of creation (annihilation) of a polariton on the lower branch with 
the wave vector k. The appearance of abnormal averages indicates a breaking of the selection 



rules connected with the gauge invariance of the system [13, H, 15]. In this problem the breaking 
of gauge symmetry arises as a result of the noninvariant initial condition due to the assumption 
that there is a condensate in the system at the initial moment t = to- 

An earlier attempt to obtain kinetic equations for polaritons excited in semiconductors by the 
external classic field was made in [|16|, 17|. However, in |L6|, [17| degeneracy of two-particle states 
was not taken into account and the abnormal distribution functions were not introduced. So 
one should expect the equations obtained in Jl6|, [HJ to possess unphysical singularity. 

Kinetic equations describing evolution of partially coherent polaritons that take into account 



the degeneracy were obtained in [18], [H| using the nonequilibrium Green's-function technique 
21] , presented by the authors in terms of functionals. They coincide with the equations 



obtained in |22[] using the method of nonequilibrium statistical operator |[23|| , and do not possess 
unphysical singularities. 



According to [T^, the kinetics of partially coherent polaritons is described in the Born 
approximation by the closed set of nonlinear integrodifferential equations for the coherent part of 
polariton field \I/ko an d the normal 

™k = - ^kol^kol 2 



and abnormal 



/k = Fk - <*k,ko*ko 



distribution functions. In the absence of quantum fluctuations, described by the functions 
and /k, the equations for them become identities, and the equation for \&ka i a is the number 
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of polariton branch) coincides with that obtained in [23| for the system of interacting coherent 
excitons and photons. In another particular case, when \l/k = and /k = 0, the equations 



obtained in [19, |22[ are reduced to the usual kinetic equation for the distribution function iVk 
(see, e.g., © [§). 

The right-hand sides of kinetic equations obtained in [|H| include terms linear in the 
constant of exciton-exciton interaction v > and the ones ~ v 2 . Terms ~ u correspond to 
the self-consistent field approximation, which neglects the higher-order correlation functions. In 
this approximation the processes of two-particle excitations of polaritons from the condensate, 
backward processes, and transformation of a created pair of noncondensate polaritons into another 
pair with the same value of total momentum are taken into account. It describes the fastest 
processes in the system and is sufficient for study of the early (before-kinetic) stage of the system's 
evolution. 

Terms ~z/ 2 take into account scattering processes in which only one polariton belongs to the 
condensate. Therefore, they differ from zero only if noncondensate polaritons exist in the system. 
Terms ~ v 2 in comparison with terms ~ v describe slower changes of the system characteristics 
and are significant only at the kinetic stage of evolution. 

Our study of time evolution of the nonequilibrium polariton system is based on the self- 
consistent field approximation. This approximation has shown a good performance in the theo- 
retical study of states of the electron-hole subsystem in semiconductors that appear just after the 
transmission of the front of the ultrashort laser pulse [^6], ^7]. The part of evolution equations 
obtained in our work for the stationary case are in many respects similar to the equations given 
in p6|, ^ p7| . Nevertheless, that slight distinctions are important and lead to significant physical 
consequencies. 

The stationary equations obtained in |26|, ||, ^] possess nontrivial solutions describing states 
of the electron-hole subsystem. The decay of these states occurs only due to incoherent relaxation 
processes. In [^, 0, |27j it is assumed that these processes do not have time to happen during the 
period of pulse action. This assumption justifies the use of the self-consistent field approximation. 
The attempt to find the steady-state solution of similar equations in the framework of the physical 
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problem of this paper leads to instabilities. Their physical nature has been discussed above. As 
noted, the appearance of the instability points leads to essential non-stationarity of the processes 
in the system and requires the return to the starting (nonstationary) equations. 

Moreover, the excitation of the great number of noncondensate modes can lead to significant 
depletion of the condensate. So to take into account this phenomenon it is necessary, along with 
equations similar to equations in [^, [27J to consider an additional equation for the condensate 
wave function ^k . Note that in p6[ 0, backward influence of the electron- hole subsystem on 
the laser radiation field was not taken into account. 



Note also |28], |29], |3(J, [y], [32] where the self-consistent field approximation allowed us to take 



into account the biexciton complex structure. 

2 System Hamiltonian 

We start from the Hamiltonian 



H = H + H- mt , 



^0 = E E Wcnoa^fcu^kaa , h k = 0^7^) + (—^) + % r 2 > ( 2 ) 

k «i,« 2 =l,2 \ Z / V Z J 



#int - 7^7 E E z/ k 1 -k 4 ^ki+k 2 ,k 3 +k 4 

V ki,...,k 4 ai,...,a 4 =l,2 

X (—iP) i^T 1 ) 0kiaift£a 2 ^«S0k4«4 • (3) 



In (|l])-(|D the following notations are used. The operators <pk a (<p£a) °^ excr ton (when a = 1) 
or photon (when a = 2) annihilation (creation) in the state with wave vector k obey Bose-type 
commutation relations: 



2 a 2 



^ki,k 2 ^ai,a 2 ) [^kiai,^k 



2 a. 2 J 



2"2 



0. 



The Pauli matrices T\, t 2 , and T3 are taken in the standard representation [20|. The frequencies 
and uj£ are given by the expressions co>k = c |k| e B 1//2 , cj^ = cu- 1 + hk 2 /2m , where c is the vacuum 
velocity of light, Ziu;- 1 " is the exciton formation energy in the band tku-^, and m is its effective mass. 



The background dielectric function e# includes the contributions from all excitations in a crystal 
except the excitons of the isolated band hu^. In the vicinity of the exciton resonance m u k 
weak frequency dependence of es can be neglected. 

The Hamiltonian ([!]) describes the system of interacting dipole-active excitons and photons 
with transverse polarization in an infinite crystal (V — > oo, V is the quantization volume) in the 
vicinity of an isolated exciton resonance m In other words, we assume that the relations 



h 



<C hujfr , AS are imposed on the photon frequency htu^, the energy of the exciton hu>^, 
and the minimum energy gap AS between the exciton band hu k and any other one. Besides that 
we suppose that the constants of exciton-photon (% = t]^) and exciton-exciton (z^ = z/_k) 
interactions are small enough: < rj-^/H <C u^r , < {p\s_/h)n <C oj£ (n is an average density 
of excitons in the system). These assumptions allow us to retain in (|2|) only resonant terms and 
treat the exciton-exciton interaction as a small perturbation. 

The quadratic part (0) of Hamiltonian ([]]) can be reduced to the diagonal form 

# = E E ^kA + Q $ ka (4) 

k a = 1,2 

by transition to the polariton Bose operators $k a and <& ka using the unitary transformation 

<£kai = E (^k)aia 2 $ka 2 > ^k = , * ^ ■ 

a 2 =l,2 yl+-^ k 

2 



L k is a function determined by the equation r/^L^ + h Lo^ — u^j L k — % = 0. In the polariton 
representation we have 

#int = 7^7 E E ^kj-ki ^ki+k 2 ,k :i +k 4 

ZK ki,... )k4 ai,...,a4=l,2 
X (^ kl k 3 ) aia3 (^k 2 k 4 ) a2a4 *k 1 a 1 $ k 2 a 2 $ k3a3^k 4 a4 > ( 5 ) 



where 

^kjk 2 = v£ ki = t/ k + 

Hamiltonian (|I|) is invariant with respect to the gauge transformation R = exp (i^N^j where 
7 is an arbitrary real parameter, N = Z)kSa = i,2 < ^ka < ^ ) ka is the operator of the total number of 
polaritons in the system. 



3 Self-consistent field approximation 

In the Heisenberg picture operators $ko(^) are governed by the equation of motion 

lhj t K«{t)= [$ka(t),#] (6) 

and average value of a dynamic quantity A is given by (A) t = Tr pA(t). Here p is the density 
matrix, which describes the polariton distribution at the initial moment of time t — to. 

The presence of the condensate in the system implies that the coherent part of the polariton 
field 

= ( $ ka) t (7) 

is nonzero. We write down the coherent part explicitly 

$ka(*)=*ka(*)+Xka(*)- (8) 

According to definition (|7|), we have (xka)t = 0- 

Using (H) we present the Hamiltonian H = Hq + H- m ^ in the form H = H\(t) + H2(t), where 
the operator H\{t) includes only linear and quadratic terms with respect to XkoSt) an d X^ a if)- 
The operator &2(t) includes products of three and four operators X'kaif) an d XyL(t)- 

Further we shall treat H 2 (t) in the self-consistent field approximation. For this purpose we 
make the formal substitution 



Xk iai (*)x£oa (0Xk 3 a 3 (*)Xk 4 a 4 (*) > 

X^ iai (^Xk 2 a 2 (*)(Xk 3 a 3 Xk 4 a 4 )t + Xk 3 a 3 (0Xk 4 a 4 (*) (Xk ia!1 Xk202^ ^kiai (*)Xk 3 a 3 (*) (Xk 2tt2 Xk 4 a 4 )i 
+Xk 1 a 1 ( t )Xk 4 a4( t )(Xk 2a2 Xk 3 a 3 )t + Xk 2 a 2 (*)Xk 3 a 3 (*) (Xk^XktcJi + Xk 2 a 2 (*)Xk 4 a 4 (*) (x£ai *k 3 a 3 ) 



t , 



Xl^aaWXksasWXkia*^) ► 

Xk 2 a 2 ( t )(Xk 3 a 3 Xk 4a4 )« + Xk 3C , 3 (*)<Xk 2 a 2 Xk4 a 4 /i ~r Xk 4 a 4 

As a result, i^ 2 (^) (and H) has the same operator structure as H\(t). Performing transformation 
from the operators Xka(t) back to the operators $ka(^) we obtain 
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ff = ^) + EE^kA + a (^k a (i) + w E E <W 2 ,k 3+k4 

k a=l,2 ZV ki,...,k 4 ai,...,a 4 =l,2 



'3/0203 



X {[^kx-k 4 (^lciks)aicx8 (^k 2 k 4 ) Q2Q4 + ^k 4 -k 3 {'Pk 1 k*) aiai (^k 2 k 
X *k iai W^k30 3 (t)($ k + 2Q2 $k 4 o 4 ) t + ^k 1 a 1 (*)*k 2 a 2 (*)( $ k 3 03 $ k404) t 

- 2$+ ai (t)n 2 a 2 m^«M^«M] + H.C.} , (9) 

where E (t) is a c- number function. 

With the help of Hamiltonian (|9]) and equations of motion (|6|) we obtain the set of equations, 
which describes time evolution of the system of partially coherent polaritons: 



*kiai(*)— y E E <Wk 2 ,k 3 +k 4 

k 2 ,k 3 ,k 4 o 2 ,a 3 ,o 4 =l,2 



X 



^ki-k 4 (^kiksJaias V^^-zkJa-icu ^ki-ks {P'k 1 'kJ a - lctA (T^k 



3/0203 



2*k 2 o 2 W^k 3 o 3 (t)^k 4 o 4 (t) + n(k 2 , a 2 ; k 3 , a 3 |t)% 4a4 (t) + -^4 a2 (t)/(k 3 , a 3 ; k 4 , a A \t) 



(10) 



+M^k iai -^k 2 o 2 ' 



n(ki,ai;k 2 ,a 2 |t) = -— E E ^k 1+ k 3 ,k 4 +k 5 

k 3 ,k 4 ,k 5 o 3 ,o 4 ,o 5 =l,2 



^ki-ks (^k^kiJa-a,, C^kgkJ aRCa + ^ki-^ (^kskJa.Q,, ("^^k 



3 /a 4 a 3 



x\n*(k 2 , a 2 ; k 5 , a 5 |t) n *(k 3 , a 3 ; k 4 , a 4 |i) + *ksaa(*)*k 4 a 4 (*) 



+-/(k 2 , a 2 ; k 3 , a 3 |t) [/ *(k 5 , a 5 ; k 4 , a 4 |t) + #k 5 o 5 (^k 4 o 4 (*) 
- {Idem [(kiai) (k 2 a 2 )]}* ; 



(11) 



^ -^(fik iai +^k 2 o 2 ; 



/(k!,ai;k 2 ,a 2 |t) = — E E 5 k 1+ k 2 ,k 3 +k 4 

K k 3 ,k 4 o 3 ,o 4 =l,2 



X 



z/ k i -k 4 (^kiks)aiaa (^£2k 4 ) Q!2 Q ! 4 + ^ki-ks (7\ik 4 ) ai a4 (^ , k 2 k 



'3/0203 



x [/(k 3 ,a 3 ;k 4 ,a 4 |t) + ^ k3a3 (t)^ k4a4 (t)] + — E E ^k 1+ k 3 ,k 4 +k 5 

k 3 ,k 4 ,k 5 o 3 ,o 4 ,Q5=l,2 



^k!-k 5 (^kik4)a,a* (^k 3 k 5 ) aaaB + ^k 4 -k 4 CPkiksLo, (^k 3 k 4 ) 



'o 3 ,o 4 



x<^ /(k2,a 2 ;k5,a 5 |t) n(k 3 , a 3 ; k 4 , a A \t) + *k, a ,(*)*k 4 a 4 (*) 



+ in*(k 2 ,a 2 ;k 3 ,a 3 |t) [/(k 5 ,a 5 ;k4,a 4 |t) + % 5Q5 (t)^ 4a4 (t)]| + Idem[(k 1 a 1 ) ^ (k 2 a 2 )] . (12) 
Here 

7i(k 1 ,ai;k 2 ,a 2 |*) = 7i*(k 2 ,a 2 ;ki,Oi|<) = ($k 1Q1 $k 2C * 2 )t - ^k^W^k.^) , 

(13) 

/(ki,ai;k 2 ,a 2 |*) = /(k 2 , a 2 ; k l5 «i|t) = ($k 1Q1 $k 2 a 2 )t - *k iai (*)*kaaa(*) ■ 
Equations (|T0|)-([T2D can be reduced to a simpler form in the physical situation when an external 
classical source is switched on and the system is prepared so that at t — to only one macroscopically 
occupied coherent mode k = k on the lower polariton branch is excited: 

*ka(*o) oc WSa^&kfo , n(ki,ai;k 2 ,a 2 |to) = /(ki, «i; k 2 , a 2 |£ ) = . (14) 

From (p!0|)-([T2]) we can find the variation of the functions ^in a (t) and (|l^) at the time moment 
t = to + dt : 

d^ ka oc 5 k ,k > dn(k 1 ,ai;k 2 ,a 2 ) = 0, df(k x , a X ] k 2 , a 2 ) oc 5 kl +k 2l 2k ■ 

Therefore, the coherent polariton wave of macroscopic amplitude is a source of quantum fluctua- 
tions, described by the function /(k 1; a.\\ 2k — k 1; a 2 \t). Substituting 

*ka(*) = ^k,k *k a(*) , n(ki, ai; k 2 , a 2 |t) = , 

/(ki, ai; k 2 , a 2 |t) = 5 kl+k2 , 2ko /(ki, a^; 2k - k 1; a 2 \t) 
into the right-hand sides of (|Tl^)-(|T2"D we obtain 

d^ka oc 5 k) k , dn(ki, a x ; k 2 , a 2 |t) oc 5 klika , t^(k ls au k 2 , a 2 \t) oc 5 kl+k2 , 2ko . 

This means that after the quantum fluctuations, described by abnormal distribution function 
/(ki, a±] 2k — ki, a 2 \t), fluctuations that are characterized by normal (usual) distribution function 
n(ki, aci] ki, a 2 \t) are excited. 
Substituting 

*ka(*) = ^k,k *k a(*) , n(ki, ai; k 2 , a 2 |t) = 5 kli k 2 n(k 1 , ati; k l5 a 2 |t) , 
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/(ki, a x \ k 2 , a 2 \t) = 5k x +k 2 ,2k /(ki, «i; 2k - k 1; a 2 \t) (15) 

into the right-hand sides of (|iT|)-(|i2|) we can easily see that the left-hand sides of these equations 
have the same structure. Therefore, solution to the (|T0|)-(|T2|) under the initial conditions ( |14]) has 
the form (|T5|). This result is nothing but a consequence of the momentum conservation law in 
the elementary acts of polariton scattering. Further simplification of (|T0|)-(|l2|) can be achieved 
by retaining in their right-hand sides only the resonant terms corresponding to lower polariton 
branch. Taking into account fll5"D and omitting fast-oscillating terms in (p!Q|)-(p^) we obtain 
A r \£ 2 (t) i 

ih j t ^o(t) = [^ko + ^ko(0]*ko(0 + ^pol (ko ' ko) -^r 1 + ^o(t)\n o (t) ; (16) 



<^n k (t) = fm 



n w 



Vl (k,ko)^^ + ^ k (t) 



1A 



(t) 

1)ol (k,k )^^ + ^(t) 



V 



(17) 



nn k + z> pol (k, k ) |VI/k ^ )|2 + ^ k (t) 

+ ^2k -k + ^p O l( 2k - k, k ) ^ k ^ )|2 + -F 2 k -k(*) 

(t) 

i/ ifk, k )^ e ^+^k(t) 



+ [1 + n k (t) + n 2ko - k (t)] • 



'pol 1 



r 



Here 



•^kW = ^E>l( k ' k l)/kiW = ^2k -kW , Fk(t) = yS^pol^klHxC*) = 4*W 
ki ki 



and 



i/ po l(k, k x ) 



^k-ki 



(l + L k ) (l+LL -k) (l +^ 1 )( 1 + L *.-k 1 
^0 + ^k-ki 



(19) 



po 



l(k,ki) 



1 + Ll) (l + L 2 ki 



are the polariton-polariton interaction constants. According to QI^D, the following initial conditions 
should be imposed on the solution to (TO)-(IH): 



"k(*o) = /k(*o) = , * ko (t )^0. 



(20) 



Further we shall assume £ = 0. 
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In the equilibrium problems .F k and are called the order parameters and are determined 
from integral equations [|1^, 571]. In the nonequilibrium problem being studied the order parameters 
depend on time. Therefore, their calculation involves solution of the set of nonlinear integrodif- 
ferential equations (|IT)D-(|T£|). 

Equations (|IT)|)-(|TBD are invariant with respect to time inversion and do not change under the 
transformation t -> -t, # ko -> # ko > / k -»• /£. 

Besides additive integrals of motion displaying the conservation laws of the average values of 
the number of particles, energy and momentum of a closed system, equations fll6|) - possess 
additional integrals of motion n k (i) — n 2ko _k(£) = const and \fk(t)\ 2 — n k(t)[l + n k (£)] = const. 
Using the initial conditions (|20|) we obtain 



Wk(*) = %k o ~k(0 



l + 4|/k(t)P 



1 



(21) 



Expressions (|2~ID enable to exclude equation (|T7|) from the set of equations ([U]) - dl8|). Note that 
representation of type (pi] ) were obtained in p6|, 0. 
It is easy to see that the transformation 



/k(i) 



■«k(*) 



l*k (0| 2; 

splits out of the set the equation for the phase of the condensate wave function. As a result, the 
set of evolution equations takes the form: 



jN (t) = ^img^(t)-N G (t) 



(22) 



ft(fi k + fi 2 k -k 



2Q 



k 0/ 



^k(t)+^ 2 k -k(*)-2^ ko (t) 



+ [z> pol (k, ko) + ^ p oi(2k - k, ko) - 2^(k , kg)] 1 k( j) }l AT (t) - 2R e e ko (t)}# k (t) 

(23) 



+ 



l + 4|rf)| I 



po 



lfcfr) 1 * 1 ^ 1 .AAo(t) + g k (t) 



where 



and 



= y ^^pol^^O^kiW = ^2k -k(*) 



k! 



Ao(t) = |* ko (t)| 2 /|* ko (0)| 2 
is the relative density of polaritons in the condensate. 
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4 Numerical results 



The set of equations (|23|) still remains rather complicated even for numerical methods. That 
is why for simplicity we neglect the dispersion of the exciton-exciton interaction constant v. In the 
vicinity of exciton-photon resonance the function L k pa 1 and, consequently, z/p Q j pa 2Vp Q \ pa v/2. 

The initial conditions (P0| ) are the same for all values of the wave vector k. Therefore, as 
follows from equations (|22]) , (^), when there is no dispersion of the polariton-polariton interaction 
constants, g depends on k via the functions f2 k + ^2k -k — 20k : 



9k = g(^k + ^2k -k - 2fi ko , t) . 



Using the identical transformation 



^E^( fi k + ^2k -k-2fi ko , t) 

k 

— fdxgiVL^ + fi 2 k -k - 2fi ko > 0<K^k + ^2k -k - 2fi k() 
K k J 

= Jdxg(x, t) ^^5(fi k + fi 2 k -k - 2fi ko -x), 



k 

and moving to the dimensionless variables 



T — Xt, w^ = A- 1 (fi k + fi 2k() _ k -2fi ko ), A = |- ' (24) 



2 



we obtain the set of equations 



Aa/" (T) = 2ImG(T).jV (T), (25) 



i-^g(w,T) = [w+Af (T)-2ReG(T)}g(w,T) 



+ Jl + 4\g(w,T)\i 



l -M,(T) + G{T) 



where 



(26) 



G(T) = Jdw g(w)g(w,T) , g(w) = 7^~j^ £ 5 ( w ^ ~ w ) ■ ( 27 ) 

Further the parameters of polaritons formed by mixing of photons and 1A excitons in CdS single 
crystal are used for numerical estimation: eg = 9.3 , Tiu 1 - = 2.55 eV , rj/fi = 1.1 x 10 14 c _1 , m± = 

13 



0.89m (m is the free electron mass), and mil = 2.85m 1 55fl . The model of isotropic parabolic 
exciton band with effective mass m 



1 /3 

m 2 ,m\\) . The effective Bohr radius of the 1A exciton 



in CdS is a ex = 28 A [53]. The exciton ionization energy is I ex = 27meV ||35|| . The value of the 
exciton-exciton interaction constant is evaluated by the formula v = Vq — (267r/3)/ ex o? ex [^] and 
makes 4.3 x 1CT 33 erg ■ cm 3 . 

At the initial stage of evolution of the system when a considerable portion of polaritons is still 
in the condensate and the number of noncondensate polaritons is rather small we can omit the 
terms containing G(T) in equations (]25|) , (|26|) . The solution to the obtained set of equations has 
the form 

JV (T) = const = 1 

(the influence of a small portion of noncondensate polaritons on the condensate is not taken into 
account), and 

1 



n(w, T) 



w(w + 2) 
1 



sinh 



sin 



^-w(w + 2) 



for 



-2 < w < 



(28) 



^/w(w + 2) 



for w < —2, w > 



w(w + 2) 

Thus, at the initial stage of evolution the excitation of noncondensate polaritons takes place 
in the range —2 < w < where the distribution function n(w) has the form of a symmetrical 
bell with a maximum at w = —1. While moving away from the above-mentioned region the 
distribution function decreases with oscillations. 

The numerical solution of the equation set (p5|), (f2~6|) was made in two stages. First, the 
explicit form of the function g(w) for different values of condensate density [^^(O)! 2 /^ and its 
wave vector k magnitude was found. Further, the function obtained was used to integrate the 
evolution equations fl2"5|), (|26|). The result of the integration for |k | = 3.6 x 10 5 cm -1 are shown 
in figures 1-3. This choice of |ko| corresponds to the exciton-photon resonance region located a 
bit below the exciton band bottom. 

The decay of the nonequilibrium polariton condensate is shown in Figure 1. Before the moment 
of dimensionless time T pa 7 the condensate depletion proceeds rather slowly. Then a drastic fall 
takes place. As a result, by the moment T pa 15 only 10% of the initial polariton number remain 
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in the condensate. The further condensate decay is accompanied by oscillations. This leads to the 
partial restoration of the condensate (up to 30% at T pa 20). The oscillations do not occur if in 
(Bp) the function G(T), describing the integral influence of all pairs of noncondensate polaritons 
with the same total momentum on each individual pair, is neglected. Thus, the occurrence of 
oscillations is due to correlation of the states of individual pairs of noncondensate polaritons. By 
the moment T pa 50 the oscillations disappear and further condensate evolution has the form of 
slow monotonous decay. These results cannot be obtained by the introduction of phenomenologi- 
cal constants in the dynamic equations. 



n(w, T) x 10 5 




20 40 60 80 100 -2 -1.5 -1 -0.5 0.5 1 

T W 

Figure 1. Dependence of a relative polariton density Figure 2. The distribution function of non- 

in the condensate on dimensionless time T. condensate polaritons at T = 10, 15, 20. 

Figure 2 depicts the distribution function of noncondensate polaritons at time moments in 
the interval where the polariton condensate undergoes the most rapid and significant changes. 
The curve corresponding to the moment of time T = 10 as well as the curve described by the 
function ( p8|) has the form of a bell with the maximum at w = —1. By the moment T = 15 the 
distribution function without changing its symmetrical form increases significantly its value and 
its maximum shifts to greater energies of noncondensate polaritons. An increase of the number 
of noncondensate polaritons by more than an order of magnitude corresponds to the presence of 
the local minimum of Mq(T) function at T pa 15. When T = 20 the distortion of the distribution 
function takes place and the area under the curve diminishes. This is in accordance with the local 
maximum of the Afo(T) at T pa 20. Thus, the distribution function is "breathing" in conformity 
with the variation of condensate density. 
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n{w,T) x 1(T 5 T 




Figure 3. The distribution function of non-condensate polaritons at different moments of 
dimcnsionless time T. The function scale is shown by the marks in the left vertical axis. 

Further evolution of the distribution function is shown in Figure 3. At the moment of time T = 
30 its contour gets the oscillating character. The oscillations arise not only in the dependence of 
distribution function on parameter w but in its dependence on dimensionless time T (unfortunately 
we could not demonstrate this in the picture). This means that if, e.g., the distribution function 
at the moment T p=* 40 possesses a local maximum at w — — 1 it can possess at this point a local 
minimum at another moment of time. Later on, the frequency oscillations become denser, smaller, 
and transform into chaotic ripples on the smooth contour of the distribution function. 

From the moment T ps 30 a short-wave wing begins to arise in the distribution function curve. 
By the moment of time T = 100 the function localizes in the interval —2 < w < 1. Very important 
is the enhancement and then the partial extinction of the sharp extra peak localized within the 
region —0.22 < w < —0.08 with a maximum at w — —0.15. The time interval where the extra 
peak arises and exists is characterized by small values of condensate density and its relatively slow 
change (see Figures 1, 3). The extra peak does not change its position with time and reaches its 
maximum at T pa 80. 
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Under the change of the initial condensate density |\I/ ko (0)| 2 /^ on ly the scale of shown depen- 
dences changes due to the corresponding transformation of dimensionless arguments (24). This 
result follows from the very weak dependence of the solution of equations (|25|), (|26|) on the initial 
condensate density via the function g(w) involved in the definitionin of the G(T)-function. Ac- 
cording to figures 2, 3, g(w,T) differs from zero only at —2 < w < 1. The numerical results 
show that in that region g(w) does not practically depend on its argument due to peculiarities of 
low polariton dispersion in the vicinity of exciton-photon resonance. This enables us to replace 
the function g(w) by g(0) in the definition (|27|) . As a result we obtain G(T) « g(0) Jdw g(w,T). 
But according to ([F?]) and (|24l) g(0) does not depend on the initial condensate density and is 



determined only by the value of wave vector k under the given crystal parameters. The value 
g(0) ~ 6.8 x 1CT 6 corresponds to the value |k | = 3.6 x 10 5 cmr x . 

5 Conclusions 

Therefore, the significant condensate depletion takes place already at the before-kinetic stage 
of the system evolution. Within a short time ~ 15A -1 the condensate density diminishes 
10 times compared to the initial density. For CdS single crystals and initial condensate density 
(10 16 -T- 10 18 ) cm" 3 the time equals to (120 -=-1.2 )ps. The short time of the polariton condensate 
decay makes problematic the possibility of observing such coherent nonlinear phenomena as, e.g., 
the phenomenon of self-induced transparency and creation of solitons in exciton spectral range 
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Note the irreversible- in-time character of the described solution of the time-reversible equations 
HI- 

We would like to point out also that in the steady-state solution l^koWl 2 = const, 
n k (t) cc 6(E k + £ 2ko _ k - 2£ ko ) , / k (t) ex 5(E k + £ 2k() _ k - 2E ko ) to equations (|T§ - (0 ) has 
been obtained ( _E k is the renormalized energy of the polariton). Substitution of this solution in 
the expressions for the total average energy and total average number of particles shows that this 
solution cannot be obtained under the initial conditions (^). But the state of the polariton system 
corresponding to this solution can be arranged by means of steady laser action on a semiconductor 
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